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NOMENCLATURE 


Cp = specific heat of air 

C = normalization constant for spectrum 

D = aperture diameter 

E^^(k) = spectrum of fluctuations, Eq. (1) 

I = light intensity, w/ cm^ 

k = wavenumber of light, cm - * 

ft - scale size of turbulence 

n = exponent in spectrum 

o = initial 

r = separation distance 

T = temperature 

u = velocity in x direction 

W = wall 

y = distance normal to boundary layer 

a = extinction coefficient, cm'* 

(3 = constant Eq. (12) 

6 = boundary-layer thickness, cm 

A = rms value 

€ = dielectric constant, air 

0 = angle from optical axis 

p = mass density 

<j) = azimuthal angle 

o = volume scattering 

ft = solid angle 

°° = free stream 
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ABSTRACT 


An analysis is presented of optical resolution quality as 
affected by aircraft turbulent boundary layers. The wind-tunnel data 
of Stine and Winovich*, was analyzed to obtain the variation of boundary 
layer turbulence scale length and mass density rms fluctuations with 
Mach number. The data gave good agreement with a mass density 
fluctuation turbulence spectrum that is either isotropic or orthogonally 
anisotropic. The data did not match an isotropic turbulence velocity 
spectrum which causes an anisotropic non -orthogonal mass density 
fluctuation spectrum. The results indicate that the average mass 
density rms fluctuation is about 10% of the maximum mass density 
across the boundary layer and that the transverse turbulence scale 
size is about 10% of the boundary layer thickness. The results indicate 
that the effect of the turbulent boundary layer is large angle scattering 
which decreases contrast but not resolution. Using extinction as a 
criteria the range of acceptable aircraft operating conditions are 
given. 
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INTRODUCTION 


The question of the quality of optical imaging through aircraft 
windows with turbulent boundary layers is important for both reconnais- 
sance and earth resources technology. In general an optical image 
recorded through such a window may be degraded by both the variation 
of the mean properties of the boundary layer in the direction parallel to 
the window and the fluctuations of mass density in the boundary layer on 
the window, since the index of refraction varies proportionally to the mass 
density. In this note we consider only the latter, and use the wind tunnel 
data of Stine and Winovich* to determine bounds for the degradation for 
flight cases. We will compare the differences between this analysis and 
previous analyses. 2 

The degradation of the image quality depends on two quantities, the 

extinction number a6 which represents the scattering of light out of the 

diffraction pattern, and the ratio of the turbulence scale size to the diameter 

of the imaging optics. If the turbulent scale size is much smaller than the 

aperture diameter, the effective resolution is insensitive to a 6 , but the 

3 

image intensity at the image plane is decreased by exp (-a6). The light 
scattered by the turbulent boundary layer raises the apparent background 
intensity. 

Analysis of Wind-Tunnel Data 

The differential cross section for scattering from a unit volume of 
a random medium is given by^’ ^ 

dcr/de = (1/4) tt{A€) 2 e" Z (1 + cos 2 0) k 4 E ee [(2k sin©/2)j (1) 
where E £e is the three dimensional spectrum of fluctuations of the dielectric 
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constant such that 


w. 

/ 

_0O 


E ec (k) dk = 1 


(2) 


The cross section for light intensity scattered out of a cone of half 
angle 0 is then given by 


o (©) = 2tt j 


(dcr/dfl) 0d0 


(3) 


and the light energy remaining in the cone is determined from 


dl (0)/dy = - I (0) a (0) 


(4) 


Using a general three-dimensional spectrum 


E ee = C * < 1 + k l > 


-n 


(5) 


Eq. (3) and (4) may be integrated. Since k » 1, 0 is always small hence 
2 sin 0/2 =0, cos © ~ 1. The result for the light intensity through a 

small aperture at the focal plane after traversing a boundary layer is: 


*[■«/■.] •-§#£/ 


( 6 ) 


•p 

If £ and ( Ae) do not vary greatly through the major portion of the 
boundary layer, we may take them as constant, equal to their average 
value. Then Eq. (6) becomes: 


/n Jl (0)/l o J 


C7T 2 k Z e'H (A€) 2 6 

a^kW ) 11 " 1 


(7) 
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The intensity "on axis", which corresponds to the image intensity is 
obtained from Eq. (7) by setting k £.& « 1, yielding: 



= - C n Z k Z Z(Ae)V 2 6/2(n-l)= -q6 


(8) 


which gives the formula for the extinction number. Eq. (7) can then be 
rewritten as 


k 2 l Z 0 C = 


|in^I(0)/l o J/in^I(e)/l o j| 


l/(n- 1) 


-1 


(9) 


There exists three unknowns in Eq. (9): i, I (0)/l Q , and n, which are 
found from the experimental data. In the experiment of Ref. 1, collimated 
light of wavelength 0. 52 jimwas passed through a wind tunnel with turbulent 
boundary layers, and refocused. Various aperture stops were used and 
the total light intensity which passed through the aperture was measured 
giving I (0)/l Q . The radius of the aperture defined 0. 

If the correct value of n is chosen, a plot of jj/n {I (0)/l Q ) 

2 

vs 0 should be a straight line. For the high Mach number, high density 
experiments it was found that the best fit of the data to a straight line 
occurred when n = 2, (see Fig. 1) which corresponds to an exponential 
correlation function of the index of refraction fluctuations, with an 
integral scale equal to H , and C = it . The intersection of the straight 
line with the abscissa then gives j£(since k is known), and the intersection 


-1 

l/(n- 1) I 


rr — - — 

The analysis of Ref. 1 assumed that n = 2. 
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Fig. 1 Theoretical and experimental* dependence of light 

intensity with aperture half angle, for two theoretical 
correlation functions for dielectric constant fluctuations. 


237 



of the line with the ordinate gives JLn (o )/i 0 J which, from Eq. (8) 

2 

gives (Ae ). For the wavelength used (0. 52pm), 

(Ae) 2 /e 2 = 4 (2. 76 x 10~ 4 )2(Ap 2 )/p 2 ( 10) 

so that the rms mass density fluctuations can be obtained. 

In Ref. 1, the Mach number was varied from 0. 4 to 2. 5 and the 
free stream density from about 1 atmosphere to 0. 1 atmosphere. In 
analyzing the data according to the above scheme, the greatest optical 
effects and hence smallest spread in the data occurred for the higher Mach 
numbers. We have analyzed only the data for the turbulent boundary layers 
with natural transition. 

RESULTS 

Figure 2 shows the resulting ratio of integral turbulence scale to 
boundary-layer thickness for 38 cases in Ref. 1. Each case is represented 
by a dot. The triangles represent arithmetic averages for each Mach 
number range. As explained previously, the largest scatter occurs for 
the lowest Mach number. However, it is clear for Mach numbers greater 
than unity that //6«0. 1, and may be somewhat larger for subsonic Mach 
numbers. These small values are initially surprising since the integral 
scale d should be about half the longitudinal velocity integral scale (if the 
turbulence is isotropic). Since it is generally believed that the ratio of 
the latter to the boundary-layer thickness is about 0. 4, we would expect 
//6 to be 0. 2. However, the scale size of the turbulence is smaller near wall, 
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the method of interpretation only yields the average value in the 
ndary layer, weighted by the fluctuations. Thus, the average value 
of £/5 is not inconsistent with our knowledge of boundary layers. 

The deduced mass density fluctuations are shown in Fig. 3, where 
we have compared the deduced mass density fluctuations to two quantities: 
the mass density difference across the boundary layer p^-Poo* anc * the 
temperature ratio across the boundary layer T w /T 0O - 1. It was assumed 
that the wall temperature of the experiment was adiabatic, with a 
recovery factor of unity. Again there is a large amount of scatter 
particularly at the lower Mach numbers. The large circles and triangles 
represent averages for each Mach number range; the former for mass 
density differences, and the latter for temperature ratio. Except for low 
Mach numbers, the ratio of ratio of A p to p ^ - p w appears to be constant, 
equal to about 0. 1, while the values based on temperature ratio decrease 
monotonically with Mach number. 

The constancy of the mass density fluctuation ratio with Mach 
number can be demonstrated approximately, as follows: assuming that 
pressure fluctuations can be neglected, then locally A p/p ~ AT/T. 

Thus the average value through the boundary layer is 


Ae_ 


w 


• o 


( 11 ) 


From Ref. 5, AT/T through the boundary layer is approximately 
constant, given by 

AT < T w/ T »- 1 > 


T x /2 (T w /Too + 1) 


( 12 ) 
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where |3 is a constant. We approximate the integral of p/p^, by the 
arithmetic average at the end points, viz: 

f d (y/6) = l/2 ( — — + A = l/2 tjt — ( 1 + ) ( 13 ) 

J Poo \P eo J 1 w A <*> 

Combining Eqs. (11, 12, and 13), 

= constant = (3 (14) 

P 00 P w 

Thus, the average mass density fluctuation ratios should be insensitive 
to Mach number. The average value of [3 from Ref. 5 is about 0. 07; 
from Fig. 3 we deduce that the average value of p is about 0. 1. This 
difference is not significant in view of the approximations made in 
evaluating Eq. (11); the main point is the constancy of Ap/(p <so -p ) 
with Mach number is consistent with measurements of temperature fluctua- 
tions in boundary layers. 

2 

To compare with previous results we note that Hufnagel used 
£/ 6 = 0 . 1 which is consistent with this analysis, however, his results 
correspond to a correlation function of the form 

-3/2 

c (r) = (1 + r 2 // 2 ) (15) 

Using Eq. (15), with Eqs. (1) and (4), the calculated dependence of I (0) 
diverges greatly from the data, see Fig. 1. In addition, the density 
fluctuations were taken to vary as 

A p = 0.2 (T w -T 0O )/l/2 (T w + T w ) (16) 
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Mass density fluctuations, deduced from the data of 
Stine and Winovich. 1 
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While this expression passes through the vicinity of the subsonic data 
points of A p/(Poo -p ) , the trend with Mach number disagrees with 
the values deduced from the experimental data; as we indicated 
previously A p/ip^ - p ) is approximately constant and does not vary 
with Mach number. 

EFFECT OF ANISOTROPIC TURBULENCE 

There is no reason to believe that the spectrum is isotropic, 
hence we have also investigated several anisotropic spectra; the only 
assumption being that one principle axis of the spectrum coincides 
with the direction of propagation. In an actual boundary layer the angle 
between principle axes of the spectrum and the spatial coordinates is 
of the order of 20°, hence the latter assumption should be quite good. 

For an anisotropic orthogonal density fluctuation spectrum, 

Eq. (5) becomes: 

E « = Ci l i 2 i 3[ 1 + ? :k i 2 ^i 2 ] < 17 > 

where, for propagation in the y direction, 

k = - k sin 0 cos <b ~ k 0 cos <i 

X t -r 

k^ = k (1 - cos 0) « k 0^/2 (18) 

k z = - k sin 0 sin (j) « k 0 sin (j) 

and dS2 - d<f> sin © d©. For £ z - £ ^ £ , then for n = 2, Eq. (7) 

x y 

becomes; 
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(19) 


■ n [i (©)/l 0 ]= 1/2 (A€) 2 / y k 2 s[l+k 2 i 2 6 2 ] 


-1 


Thus the turbulence scattering angle depends on the scale size in the plane 
normal to the path hut the extinction depends on the scale size in the 

direction of propagation {£ y ). In this case it is not possible to isolate 
separately the effects of scale size i and density fluctuations since they 
appear as a product. 

Another example of an anisotropic spectrum corresponds to an 
isoenthalpy flow in the x direction. Then Cp A T = UAu, and the spectrum 
for temperature {or density) fluctuations is the same as for A u, e. g. 

E = E^. For isotropic velocity fluctuations, we can also use an 
exponential velocity fluctuation correlation function, for which^ 

E ee = E n = 2 7r- Z t 5 (k 2 -k 2 ) (1 + k 2 i 2 ) 1 (20) 


The dependence on aperture angle 0 of light intensity is again obtained by 
integration of Eq. (1) and (4): 


2 

- */ 2 (-¥-) 


? 7 7 7 

kr /5(1 + K0t; 


'[ 


1 - (2 + 2k z t z e z 


r 


( 21 ) 


This form is similar to Eq. (1), except that there is some curvature in the 
plot of l/£n|l(0)/lj vs 0 2 for k0/ <5. This curvature is not evident in 
the data; hence the assumption of isoenthalpy and isotropic velocity 
turbulence which leads to Eq. (21) is rejected. 
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A third possible anisotropic form for E corresponds to iso- 
enthalpy flow, but with different scale sizes in the 3 orthogonal directions. 
A form of the velocity spectrum which satisfies incompressible continuity 
is : 


ij 


= A 


I . I . 
i 3 


6,-Yl k. 




k. k. I. 
i J 1 


I 2 l 

J 


f(k ; 


V 


( 22 ) 


For example, if the velocity correlation function is exponential, 
in the three orthogonal directions, then 

E ee =E ll = Z *' 2 V2 i 3f k 2 , / +(k 3 , 3 ): j[ 1+ £ k ! 2 '/] 3 (23 > 

for which Eqs. (1) and (4) may be approximately integrated to obtain 


■ n I(0)/I O 1= 


- i/2(^ e 2 )k 2 i 


l i z l : 


6 1 


-4 


j^i + k 2 e 2 f 2 j ji - 


(2+ 2k 2 f 2 e 2 ) 


- 

) 


(24) 


where 


f 2 = 3/4 i 2 2 + 1/4 


(25) 


Again, Eq. (24) exhibits curvature near the origin, which is not present 
in the data. Thus we conclude that either Eq. (7) with n = 2, or its aniso- 
tropic form Eq. (19) are most consistent with the data; and that the 
scattered light is not consistent with isoenthalpy flow with either an 
isotropic or anisotropic velocity spectrum. 

APPLICATIONS TO FLIGHT 

For the usual cases that the aperture of the imaging device is 
much larger than the turbulence scale size, the primary effect on the 
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diffraction pattern from a distant point source will be attenuation of the 
light by turbulence scattering, but the resolution will not be changed 
appreciably . Typical diffraction patterns are shown in Fig. 4. For 
example the case a 6 = 0. 50 corresponds to a boundary layer 11 cm thick, 
a Mach number of 0. 8, and an altitude of 9 km. While the theoretical 
resolution is insensitive to the extinction, excessive extinction causes 
loss of contrast, since the scattered light enters adjacent resolution 
cells. Thus, we may define some approximate criteria: if aS >2, then 
the ability to image will be poor, but if a 6 <0. 4, the ability to image 
should be quite good. Figure 5 shows these approximate boundaries, for 
three boundary layer thicknesses, 1, 3, and 10 cm. 

From Fig. 5, it can be seen that the seeing effects are very 
sensitive to Mach number in the vicinity of Mach one, but insensitive to 
Mach number at high Mach numbers. This latter effect is caused by the 
fact that 1 -p w /Poo goes asymptotically to unity with increasing Mach 
number. For all Mach numbers, the "seeing" is sensitive to the boundary- 
layer thickness, hence there is always a gain in keeping the boundary layer 
over the window of imaging device as thin as possible. Other techniques 
may also be useful for improving the imaging, such as cooling the window 
to reduce the temperature excess, and hence mass density fluctuations. 

One note of caution: these results are based on a smooth window flush 
with the aircraft skin. Window moldings, recesses, etc. could degrade 
the imaging quality by increasing either the optical path through the 
turbulence layer or the scale size of the turbulence. 
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Fig. 4 Diffraction patterns. 
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Fig. 5 Imaging quality boundaries for various boundary-layer 
thicknesses, flight Mach numbers, and altitudes. 
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